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1. INTRODUCTION 

Rare earth elements [REEs] include a series of lanthanide properties. Because scandium [Sc] and yttrium [Y] exhibit 
[Ln] elements, the elements with atomic numbers 57 to 71; similar properties to the lanthanide family, they are also 

and they have similar electronic structures and chemical considered rare earth elements [1], Two to three groups of 




ABSTRACT 


Fifteen stream sediment samples were collected in Bwanebwa area, situated 
southwestern part in the central domain of the Pan-African Fold Belt in 
Cameroon. Samples were analyzed by Aqua Regia - Inductive Coupled Plasma - 
Mass Spectrometry analytical technique. The chemical analyses revealed a 
wide variability of lanthanide concentrations [IREE = 1471.56-7674.2 ppm] 
displaying a normal distribution pattern. Lanthanides show a high enrichment 
in La, Ce, Pr, Nd, and the proportion enrichment of REEs is still similar with 
increasing atomic number. The order of the average concentrations of REEs is 
Ce > La > Nd > Pr > Sm > Gd > Dy > Er > Tb > Ho >Yb > Eu > Tm > Lu, which 
was different to that found in the Oddo-Harkins rule, earth’s crust, upper 
continental crust, South China Sea, mean REE values in oceanic crust, Redang 
Island Marine Sediment in Terengganu Coastal Waters, Bayan Obo ores, Edea, 
suggesting that REEs set up a complex group. The LREE/HREE ratio [31.57- 
35.24, av. = 32.82] showed that the content of LREE [1425-7449.8, av. = 3699] 
is significantly higher than that of HREE [46.56-224.4, av. = 112.07], LREEs 
account for 97.06% of the total REEs and this percentage is in agreement with 
the trend that was observed in the earth’s crust and in the Bayan Obo ores. 
Chondrite-normalized REE distribution patterns are identical, indicating the 
consistency of geochemical distribution of REEs in stream sediments. They are 
strongly fractionated, with LREEs significantly fractionated [LaN/YbN = 11.74- 
27.93, av. = 21.04] relative to HREEs [Gd N /Yb N = 7.40-16.69, av. = 13.35], 
implying that HREEs are more soluble and more complex than middle or LREE 
and are more strongly absorbed on most substrates. Positive Ce anomaly 
[Ce/Ce* = 1.04-1.07] and negative Eu anomaly [Eu/Eu* = 0.097-0.17] were 
observed, indicating that differentiation occurred between Ce, Eu and other 
lanthanides. Therefore the positive Ce anomaly recorded represents an 
enrichment of Ce compared with its neighboring elements. A five-factor model 
of the REEs that accounts for 1.998% of the total data cumulative variance 
have been observed, with the communality values [0,996-1] reflecting 
stronger community structure. Principal component analysis results in five 
components that explain 100% of the total variance, with the placement of 
elements in three major principal components, in which elements scattering in 
the rotated space displaying their mutual association. The results provide 
evidence for similar input sources and common geochemical characteristics of 
REEs. Lanthanide concentrations in Bwanebwa area were of geogenic origins 
and not influenced by anthropogenic sources. The abundance of REEs is 
related to the lithology made up of granites and associated pegmatites. Cerium 
shows the highest average content of REEs, followed by La, Nd, Pr, Gd, Sm, 
hence pattern of Ce distribution shows higher than the other lanthanides. 
LREEs [La, Ce, Pr, Nd] and Ga fall above the upper background threshold limit 
and are considered as anomalous. Therefore the results of this work serve as 
guide for lanthanides exploration in Bwanebwa area. 
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REEs are generally distinguished into light rare earth 
elements (LREEs], heavy rare earth elements (HREEs], and 
sometimes middle rare earth elements (MREEs] based on 
their atomic mass and effective ion radius [2] and [3], The 
LREEs or elements of the cerium subgroup are the lower 
atomic weight elements, (atomic mass less than 153 g/mol 
and ionic radius > 95 pm]. This subgroup includes elements 
from lanthanum to europium. The HREEs or elements of the 
yttrium subgroup with atomic mass above 153 g/mol, 
encompass lanthanides from gadolinium to lutetium, 
including yttrium, which is similar in its chemical properties 
and whose mass is equal to 88.9 g/mol [4], Sm to Ho are the 
elements of the MREEs. REEs generally occur as trivalent 
ions, Ce can be also in tetravalent state, and Eu is in divalent 
state. 

Lanthanides are highly electropositive and are 
predominantly trivalent (Ln 3+ ), with the exception of cerium 
(Ce 4+ ) and europium (Eu 2+ ) in some environments. The entry 
of new electrons in the 4f orbitals, when the atoms have fully 
occupied the 6s orbital increases the electrostatic attraction 
between the N shell and the nucleus. This leads to a 
reduction of the REE 3+ ionic radius with increasing Z, a 
phenomenon which is known as lanthanide contraction [5] 
and [6], Since most Ln possess similar atomic radii and 
oxidation states, they can substitute for each in various 
crystal lattices. This capability of substitution leads to 
multiple REE occurrences within a single mineral and has 
resulted in a wide distribution within the Earth’s crust [7], 
REEs are found in a wide range of mineral types, including 
halides, carbonates, oxides, phosphates and silicates. The 
abundance of lanthanides within Earth’s crust varies widely 
across individual REEs, ranging from the most abundant at 
66 ppm of cerium (exceeding other important metals 
including copper: 27 ppm and lead: 11 ppm) to 0.28 ppm for 
thulium [8] and [9], 

During the last few decades, REEs have become important 
geochemical tracers in order to understand and describe the 
chemical evolution of the earth’s continental crust [10]-[13], 
Rare earth elements have been used, as analogues for 
actinide elements, in studies related to radioactive waste 
disposal in order to demonstrate their general immobility in 
weathering environments [14], The comparison among REE 
facilitates the normalization of analyses to reference 
standards such as chondrite [15], Technological 
developments have increased the use of REE in the defense, 
aerospace, medical and automotive industries. Lanthanides 
are essential components of catalysts, high-strength 
magnets, super-alloys, display technology and lasers [10], 
[16], [17], Their concentration has been proven to occur in a 
wide range of geological settings [18], [19], 

World resources of lanthanides are contained primarily in 
the minerals bastnasite, monazite, loparite, and in ion- 
adsorption clays [20], Bastnasite deposits in alkaline rocks 
and carbonatites of China and the United States constitute 
the largest percentage of the world's economic resources, 
while monazite deposits in Australia, Brazil, China, India, 
Malaysia, South Africa, Sri Lanka, Thailand, and the United 
States make up the second largest segment [21], Brazilian 
deposits account for about 1% of the world reserve [22], 
Currently, the use of REEs is closely associated with high- 


tech industry. The increasing use of REEs in industrial 
processes suggests that REEs release into the environment is 
likely to increase in the future with potential impacts on 
human health [23], 

Despite vast global distribution, REEs are mainly mined in 
China. In 2011, China produced over 90% of the world’s rare 
earth supply while owning only 23% of the world’s total 
reserves [24], Due to the increase of rare earth mineral 
mining, China’s rare earth reserve has drastically decreased 
from 75% in 1970 to 23% in 2011 [25] and [24], The surface 
mining and heap leaching of China’s unique ion-adsorption 
rare earth resources have caused severe environmental 
damage, such as soil erosion, pollution, and acidification. 
Moreover, REEs have long been used in China as additives in 
fertilizers and as growth promoters in livestock feed. The 
long-term effects and damage of such uncontrolled discharge 
of REE-based chemicals into the environment remain to be 
determined. Lanthanides have been characterized neither as 
essential elements for life nor as strongly toxic elements in 
the environment [26], Although the environmental toxicity 
of REEs is largely unknown, environmental contamination 
has already been found in some mineralized areas as well as 
soils that are affected by the long-term application of sludge. 
Several negative effects of REEs on organisms have been 
reported. 

Little are known about lanthanides in Cameroon, meanwhile 
preliminary studies have been done by [27] who reported 
for the first time the distribution of rare earth elements from 
the metamorphic rocks of the Paleoproterozoic Nyong unit 
situated at the west of Congo Craton in Cameroon. The 
present study is carry out in Bwanebwa area situated around 
the central domain of the Pan-African Fold Belt in Cameroon. 
The aim of this work is to study the chemical distribution 
and fractionation of lanthanides in stream sediments 
samples collected in Bwanebwa area. 

2. Regional Geology 

The Bwanebwa region is located in the southwestern part of 
the central domain of the Pan-African Fold Belt (PAFB) in 
Cameroon (Fig.l). The PAFB is a major Neoproterozoic 
Orogen linked to the Trans-Saharan Belt of western Africa 
and to the Brasiliano Orogen of NE Brazil [28] and [29], In 
Cameroon the PAFB is made up ofthree main domains [30]- 
[36], (Fig. 1): (1) the southern domain, which corresponds to 
the northern edge of the Congo Craton, formed by 
Neoproterozoic metasediments deposited in a passive 
margin environment and metamorphosed under high 
pressure conditions at 616 Ma; (2) the northern domain 
consisting of subordinate 830 Ma old metavolcanic rocks of 
tholeiitic and alkaline affinities; (3) the central domain that 
contains the present study area, is situated between the 
Sanaga shear zone (SSZ) and the Betare-Oya shear zone 
(BOSZ) [37] to the south, and the Tibati-Banyo fault (TBSZ) 
to the north [38], This domain is dominated by a NE-SW 
elongated regional-scale plutonic complex intrusive into a 
Paleoproterozoic basement and locally covered by 
Cretaceous deposits (Mbere and Djerenr basins) and by 
Cenozoic volcanic rocks of the Cameroon Volcanic Line 
(CVL). The study area is made up of granites and the 
associated pegmatites. 
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Figurel. Geologic map of Cameroon, modified from [31], 
[35], [39]-[41], showing major lithotectonic units and shear 
zones: [1] southern domain; [2] central domain; [3] 
northern domain; CCSZ, Central Cameroon shear zone; TBSZ, 
Tchollire-Banyo shear zone; SSZ, Sanaga shear zone; BOSZ, 
Betare- Oya shear zone. The location of the study area is 
marked by a red box. 


3. Methods of investigation 

Fifteen stream sediment samples were collected in the study 
area. Prior to the sampling process, the areas were surveyed, 
in order to see the favorable sites for sampling. Sampling 
points are found along the stream banks where the flow 
velocity is low and favors deposition [area with no limited 
sediment accumulation]. The GPS coordinates of the chosen 
sites were then taken and recorded on the field notebook. 
Stream sediment samples were collected downstream at 
interval of about 300m, at depth of about 30-50cm to avoid 
high contents of Fe-Mn oxides coating. During the panning 
process, large particles were removed with the hands. 
Panning was carefully done for at least 30-45 minutes and 
the obtained concentrates collected. The collected 
concentrates were put into polythene bags and labelled. The 
samples were then air dried, weighed on clean paper sheets 
to avoid contamination at the geological laboratory of the 
University of Buea for a period of one week. After drying, 
about 5g of the bulk samples were packaged and sent to 
Activation Laboratories [ACTLABS] in Canada for chemical 
analysis. Chemical analysis was performed on the fifteen 
stream sediments samples using Aqua Regia - Inductive 


Coupled Plasma- Mass Spectrometry [AR-ICP-MS] analytical 
technique. 

4. Results 

4.1. Geochemical data 

The chemical composition of lanthanides of stream 
sediments from Bwanebwa area is represented in table 1. 
The chemical analyses revealed a wide variability of REE 
contents with IREE in the range of 1471.56 to 7674.2 ppm 
with an average of 3811 ppm (Table 1], The overall stream 
sediment samples show a high enrichment in La [332-1750 
ppm, av. = 863.8 ppm], Ce [703-3758 ppm, av. = 1846.93 
ppm], Pr (76.11-392 ppm, av. = 196.89 ppm] and Nd (269- 
1440 ppm, av. = 715.53 ppm]. Only samples EB05 and EB07 
show relative enrichment in Yb (4.6 and 4 ppm respectively; 
Table 1], The concentrations of individual REEs tend to 
decrease with increasing atomic number, REEs with even 
atomic numbers are more frequent than their neighbors 
with odd atomic numbers. The order of the average 
concentrations of REEs in stream sediments around 
Bwanebwa area is as follows (in decreasing order]: 
Ce>La>Nd>Pr>Sm>Gd> Dy>Er>Tb>Ho>Yb>Eu>Tm>Lu. 

The LREE/HREE ratio ranged largely from 31.57 to 35.24 
with an average of 32.82. This ratio showed that the content 
of LREE (1425-7449.8, av. = 3699] is significantly higher 
than that of HREE [Gd-Lu] that yields between 46.56 and 
224.4 with 112.07 as average. The content of LREEs [La-Eu] 
accounted for 97.06% of the total REEs content in the 
investigated stream sediments, just as they are in the earth’s 
crust [3], This percentage is in agreement with the trend 
observed in the Bayan Obo ores [42], 

The distribution of REEs in stream sediments normalized to 
upper continental crust [UCC] after [9] are illustrated in Fig. 
2. The REEs normalized patterns are strongly fractionated 
(LaN/YbN = 11.74-27.93], They show LREEs significant 
fractionation (LaN/YbN = 11.74-27.93, av. = 21.04; CeN/SniN = 
1.12-2.48] relative to HREEs (Gd N /Yb N = 7.40-16.69, av. = 
13.35], They display a pronounced negative Eu anomaly 
(Eu/Eu* = 0.097-0.17], a very slightly positive Ce anomaly 
(Ce/Ce* = 1.04-1.07] and a very weak positive Nd anomaly 
(Nd/Nd* = 0.97-1.46], Also, REEs distribution patterns are 
similar and parallel. 



Figure2. Upper continental Crust normalized REE patterns. 
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TABLE I Lanthanides Characteristics of Stream Sediments from Bwanebwa 


REE 

EB01 

EB02 

EB03 

EB04 

EB05 

EB06 

EB07 

EB08 

EB09 

EB10 

EB11 

EB12 

EB13 

EB14 

EB15 

Mean 

La fppm) 

333 

835 

710 

1159 

1753 

1522 

1350 

940 

455 

405 

425 

505 

440 

940 

1185 

863.8 

Ce 

703 

1792 

1503 

2492 

3758 

3290 

2879 

2010 

976 

855 

927 

1062 

945 

2020 

2492 

1846.9 

Pr 

76.11 

191 

159 

265 

392 

352 

308 

218 

107 

91.3 

100 

115 

100 

215 

264 

196.89 

Nd 

270 

704 

579 

967 

1441 

1270 

1121 

794 

385 

325 

362 

420 

370 

775 

950 

715.53 

Sm 

41.79 

82 

66 

100 

101 

102 

100 

91.9 

44.5 

48 

47.8 

49 

44.3 

81 

101 

73.35 

Eu 

1.1 

2.4 

1.7 

3.2 

4.8 

4.5 

4 

2.7 

1.2 

1 

1 

1.3 

1.4 

2.6 

3.4 

2.42 

Gd 

27.1 

66.8 

54.73 

94 

138 

124 

105 

76.8 

35.8 

34.2 

37.3 

38.4 

36 

73 

92 

68.88 

Tb 

2.41 

5.7 

4.7 

8.22 

12.4 

10.9 

9.3 

6.6 

3.1 

3.2 

3.4 

3.3 

3.1 

6.3 

8.1 

6.05 

Dy 

10.12 

23.1 

17.5 

33 

49 

44.3 

38 

26.2 

12 

12.7 

13.6 

13.2 

13 

25 

32 

24.18 

Ho 

1.5 

3 

2.2 

4.12 

6.3 

5.9 

5 

3.4 

1.6 

1.6 

1.8 

1.7 

1.8 

3.1 

4.1 

3.14 

Er 

3.1 

6 

4.3 

8.23 

12.5 

11.4 

10 

9 

3.4 

3.2 

3.5 

3.6 

3.7 

6.4 

8 

6.42 

Tm 

0.3 

0.51 

0.4 

0.7 

1.1 

1 

0.9 

0.6 

0.3 

0.3 

0.3 

0.3 

0.4 

0.5 

0.7 

0.55 

Yb 

1.83 

2.5 

1.64 

3 

4.6 

4.8 

4 

3 

1.4 

1.5 

1.4 

1.7 

2 

2.3 

3.1 

2.56 

Lu 

0.2 

0.3 

0.2 

0.3 

0.5 

0.5 

0.4 

0.3 

0.1 

0.2 

0.1 

0.2 

0.2 

0.2 

0.3 

0.27 

EREE 

1471.56 

3714.31 

3104.37 

5137.77 

7674.2 

6743.3 

5934.6 

4182.5 

2026.4 

1782.2 

1924.2 

2214.7 

1960.9 

4150.4 

5143.7 

3811 

LREE 

1425 

3606.4 

3018.7 

4986.2 

7449.8 

6540.5 

5762 

4056.6 

1968.7 

1725.3 

1862.8 

2152.3 

1900.7 

4033.6 

4995.4 

3699 

HREE 

46.56 

107.91 

85.67 

151.57 

224.4 

202.8 

172.6 

125.9 

57.7 

56.9 

61.4 

62.4 

60.2 

116.8 

148.3 

112.07 

LREE/ 

HREE 

30.61 

33.42 

35.24 

32.90 

33.20 

32.25 

33.38 

32.22 

34.12 

30.32 

30.34 

34.49 

31.57 

34.53 

33.68 

32.82 

fLa/Yb)N 

11.74 

21.55 

27.93 

24.92 

24.59 

20.46 

21.77 

20.22 

20.97 

17.42 

19.59 

19.17 

14.19 

26.37 

24.66 

21.04 

fGd/Yb]N 

7.40 

13.36 

16.69 

15.67 

15 

12.92 

13.12 

12.8 

12.79 

11.4 

13.32 

11.29 

9 

15.87 

14.84 

13.35 

fCe/SmlN 

1.12 

1.46 

1.52 

1.66 

2.48 

2.15 

1.92 

1.46 

1.46 

1.19 

1.29 

1.44 

1.42 

1.66 

1.64 

1.59 

fLa/SmlN 

1.08 

1.38 

1.46 

1.57 

2.35 

2.02 

1.83 

1.39 

1.39 

1.14 

1.20 

1.40 

1.35 

1.57 

1.59 

2.18 

Eu/Eu* 

0.13 

0.13 

0.12 

0.14 

0.17 

0.16 

0.16 

0.13 

0.12 

0.10 

0.10 

0.12 

0.14 

0.14 

0.14 

0.15 

Ce/Ce* 

1.04 

1.06 

1.05 

1.06 

1.07 

1.06 

1.05 

1.05 

1.04 

1.05 

1.06 

1.04 

1.06 

1.06 

1.05 

1.43 

Nd/Nd* 

0.97 

1.14 

1.14 

1.20 

1.46 

1.36 

1.29 

1.13 

1.13 

0.99 

1.06 

1.13 

1.12 

1.19 

1.18 

1.6 


TABLE II Description Statistics of Rees Concentration in Stream Sediment from Bwanebwa 


REEs 

Min 

Max 

Mode 

Range 

Mean 

GM 

Median 

Std Dev 

CV (%) 

Variance 

Skewness 

Kurtosis 

IQR 

Ql 

Q3 

UBT 

La 

333 

1753 

940 

1420 

864 

755.85 

835 

451 

52.21 

203354 

0.57 

-0.78 

745 

440 

1185 

345.93 

Ce 

703 

3758 

2492 

3055 

1847 

1614.44 

1792 

969 

52.45 

938341 

0.59 

-0.75 

1547 

945 

2492 

524.79 

Pr 

76.1 

392 

100 

315.89 

196.9 

172.68 

191 

101.8 

51.72 

10368.4 

0.55 

-0.86 

165 

100 

265 

172.44 

Nd 

270 

1441 

- 

1171 

715.5 

626.47 

704 

371.9 

51.98 

138343.1 

0.56 

-0.79 

597 

370 

967 

336.5 

Sm 

41.79 

102 

100/101 

60.21 

73.35 

69.07 

81 

25.16 

34.3 

633.03 

-0.09 

-1.97 

52.2 

47.8 

100 

144.83 

Eu 

1 

4.8 

1 

3.8 

2.42 

2.09 

2.4 

1.316 

54.38 

1.732 

0.55 

-1.03 

2.2 

1.2 

3.4 

22.46 

Gd 

27.1 

138 

- 

110.9 

68.88 

60.59 

66.8 

35.48 

51.51 

1258.59 

0.6 

-0.77 

58 

36 

94 

99.96 

Tb 

2.41 

12.4 

3.1 

9.99 

6.049 

5.32 

5.7 

3.151 

52.09 

9.928 

0.68 

-0.61 

5.02 

3.2 

8.22 

30.2 

Dy 

10.12 

49 

- 

38.88 

24.18 

21.28 

23.1 

12.64 

52.27 

159.76 

0.69 

-0.7 

20 

13 

33 

58.1 

Ho 

1.5 

6.3 

1.6/1.8 

4.8 

3.141 

2.78 

3 

1.626 

51.76 

2.644 

0.77 

-0.6 

2.42 

1.7 

4.12 

25.62 

Er 

3.1 

12.5 

- 

9.4 

6,422 

5.70 

6 

3.23 

50.3 

10.434 

0.6 

-1.02 

5.5 

3.5 

9 

29.07 

Tm 

0.3 

1.1 

0.3 

0.8 

0.554 

0.50 

0.5 

0.2721 

49.11 

0.074 

0.86 

-0.45 

0.4 

0.3 

0.7 

5.48 

Yb 

1.4 

4.8 

1.4/3 

3.4 

2.585 

2.37 

2.3 

1.141 

44.16 

1.303 

0.85 

-0.41 

1.46 

1.64 

3.1 

19.92 

Lu 

0.1 

0.5 

0.2 

0.4 

0.2667 

0.24 

0.2 

0.1234 

46.29 

0.0152 

0.74 

-0.03 

0.1 

0.2 

0.3 

66.14 


GM: Geometric mean; Std Dev: Standard deviation; CV: Coefficient of variation or relative standard deviation 


RSD); 


IQR: Interquartile range; 


Ql: First quartile; Q3: Third quartile; UGT: Upper Background Threshold 


WwtWt 
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Figure3. Cumulative Frequency Plots of Rare Earth Elements 


4.2. Data distribution patterns 

Table 2 presents the descriptive statistics of lanthanides in stream sediments. A summary of the main statistical parameters 
(mean, range, median, and standard deviation] was also shown. The geometric mean (GM) were used to describe the average 
concentrations [43], The stream sediments content the LREEs (La, Ce, Pr, Nd, Sm, Eu] and the HREEs (Gd, Tb, Dy, Er, Y, Lu). 
They are characterized by high abundance of Ce (703-3758 ppm], La (333-1753 ppm], Nd (269-1440 ppm], relative abundance 
of Pr (76.11-392 ppm], and depletion of HREEs. The average values (GM] are closer to the median values for all the REEs 
indicating a symmetric data distribution pattern RE E (Table 2], Moreover, the very low positive skewness values suggest that 
individual REE datasets were very weakly right-skewed model. Also, the very low negative kurtosis values indicate a very 
slightly light-tailed distribution. However, the studied REEs still follow a normal distribution. The graphical inspection of 
original data distribution for all REEs shows a clear normal distribution (Fig. 3], The standard deviation values were high for 
some elements, particularly for LREEs (La to Nd], and Gd (Table 2], The coefficient of variation (CV] or the relative standard 
deviation (RSD] can be used to compare in relative terms the variability of the same property under similar values of variance 
and different means. Low CV values correspond to a spatially homogeneous distribution of REE concentrations, whereas high 
CV values indicate a non-homogenous surface distribution [44], The relative standard deviation values are moderate for the 
studied lanthanides, greater to 50% (Table 2], indicating a moderate variability. The spatial variability of Sm, Tm, Yb and Lu, 
however, was low with the CVs only fluctuating from 34.3 to 49.11%. 

4.3. Multivariate statistical analysis (MSA) 

In an attempt to identify REE associations, Pearson’s correlation matrix and factor analysis was carried out to investigate the 
inter-element relationship. The analysis factors related variables into principal associations on the basis of their mutual 
correlation coefficients [45], 

4.3.1. Correlation analysis 

The correlation analysis (Table 3] between REEs exhibits an unusually case where all the REEs display higher positive 
correlation (r = 0.9144-0.9998] with their respective correspondents. Only correlation values between Sm and Tm, Yb and Lu 
that are slightly low with r values ranged from 0.823 5 to 0.8885. These findings provide evidence for similar input sources and 
common geochemical characteristics of the elements. 

4.3.2. Factor Analysis 

Table 4 shows a five-factor model of the REEs that accounts for 1.998% of the total data cumulative variance. Factor 1 (La, Ce, 
Pr, Nd, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu] account for 0.363% of the total data variance. This elemental grouping is relatively 
and positively loaded with lithophile elements with respect to the LREE and HREE. The significant loading suggests REE 
fractionation in the study area that involved both the LREEs and HREEs. The strong loading of Ho (0.616], Er (0.609], Tm 
(0.659], Yb (0.747] and Lu (0.828] establishes their mutual relationship as HREEs. Factor 2 (La, Ce, Pr, Nd, Eu, Gd, Tb, Dy, Ho, 
Er, Tm] account for 0.336% ofthe total variance. The strong loading of La (0.637], Ce (0.64], Pr (0.63], Nd (0.632] reflects their 
association as LREEs while the strong loading of Gd (0.627], Tb (0.64] and Dy (0.623] indicates their relationship as MREEs. 
Factor 3 (La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er] account for 0.29% of the total variance, and Sm (0.805] is strongly loaded. 
Factor 4 and Factor 5 are not loaded and they account for 0.007 and 0.003% of the total data variance respectively. The 
communality values range from 0,996 to 1 (Table 4] and indicate a common variance. 


TABLE III Pearson’s Correlations between Ree Concentrations 



La 

Ce 

Pr 

Nd 

Sm 

Eu 

Gd 

Tb 

Dy 

Ho 

Er 

Tm 

Yb 

Lu 

La 

1 














Ce 

0.9998 

1 













Pr 

0.9995 

0.9997 

1 












Nd 

0.9995 

0.9998 

0.9998 

1 











Sm 

0.9347 

0.9329 

0.9366 

0.9351 

1 










Eu 

0.9934 

0.9934 

0.9943 

0.9937 

0.9344 

1 









Gd 

0.9985 

0.9990 

0.9990 

0.9989 

0.9349 

0.9929 

1 
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Tb 

0.9975 

0.9980 

0.9975 

0.9974 

0.9256 

0.9909 

0.9992 

1 







D V 

0.9957 

0.9965 

0.9964 

0.9962 

0.9241 

0.9932 

0.9983 

0.9990 

1 






Ho 

0.9910 

0.9920 

0.9924 

0.9921 

0.9120 

0.9936 

0.9945 

0.9956 

0.9982 

1 





Er 

0.9768 

0.9778 

0.9800 

0.9801 

0.9260 

0.9817 

0.9818 

0.9803 

0.9822 

0.9833 

1 




Tm 

0.9808 

0.9816 

0.9814 

0.9818 

0.8885 

0.9883 

0.9837 

0.9856 

0.9895 

0.9945 

0.9789 

1 



Yb 

0.9508 

0.9522 

0.9543 

0.9537 

0.8662 

0.9733 

0.9569 

0.9571 

0.9661 

0.9780 

0.9716 

0.9847 

1 


Lu 

0.9144 

0.9149 

0.9149 

0.9159 

0.8235 

0.9322 

0.9195 

0.9207 

0.9284 

0.9397 

0.9281 

0.9507 

0.9709 

1 


TABLE IV Varimax Rotated Factor Loading and Communalities 


Variable 

Factor 1 

Factor 2 

Factor 3 

Factor 4 

Factor 5 

Communality 

La 

0.541 

0.637 

0.547 

-0.04 

0.006 

0.999 

Ce 

0.543 

0.64 

0.542 

-0.046 

0.007 

1 

Pr 

0.544 

0.63 

0.55 

-0.054 

0.015 

0.999 

Nd 

0.545 

0.632 

0.547 

-0.055 

0.007 

0.999 

Sm 

0.413 

0.423 

0.805 

-0.042 

0.024 

1 

Eu 

0.596 

0.579 

0.547 

-0.044 

0.084 

0.998 

Gd 

0.553 

0.627 

0.544 

-0.063 

0.006 

1 

Tb 

0.559 

0.64 

0.523 

-0.065 

0.004 

0.999 

Dy 

0.58 

0.623 

0.518 

-0.065 

0.029 

0.998 

Ho 

0.616 

0.607 

0.492 

-0.073 

0.057 

0.998 

Er 

0.609 

0.533 

0.546 

-0.214 

0.017 

1 

Tm 

0.659 

0.59 

0.448 

-0.078 

0.082 

0.996 

Yb 

0.747 

0.48 

0.432 

-0.108 

0.113 

0.999 

Lu 

0.828 

0.397 

0.391 

0.004 

-0.058 

1 

Variance 

5.0881 

4.7057 

4.0601 

0.0948 

0.0352 

13.984 

% Var 

0.363 

0.336 

0.29 

0.007 

0.003 

0.999 

Cum % 

0.363 

0.699 

0,989 

0,996 

0,999 

1,998 


4.3.3. Principal Component Analysis (PCA) 

PCA (Table 5} results in five components that explain 100% of the total variance (eigenvalue] for the studied samples. The first 
component (PCI] that accounts for 100% of the total variance contains Ce (0.802], this element is strongly loaded reflecting 
high Ce enrichment relative to the other REEs. The second (PC2], third (PC3], fourth (PC4] and fifth (PC5] components 
accounting for 0% of the total variance respectively. The component plot (Fig. 4] obtained after component analyses show the 
placement of elements in three major principal components (PCI, PC2, PC3], in which elements scattering in the rotated space 
displaying their mutual association. 


TABLE V Rotated Component Matrix 


Variable 

PCI 

PC2 

PC3 

PC4 

PC5 

La 

0.397 

0.074 

-0.726 

-0.536 

- 0.111 

Ce 

0.852 

0.199 

0.129 

0.446 

0.112 

Pr 

0.09 

-0.147 

0.082 

0.011 

-0.798 

Nd 

0.327 

-0.513 

0.54 

-0.56 

0.105 

Sm 

0.021 

-0.818 

-0.39 

0.405 

0.109 

Eu 

0.001 

-0.003 

0 

-0.004 

-0.041 

Gd 

0.031 

-0.031 

0.057 

0.175 

-0.387 

Tb 

0.003 

0.006 

0.004 

0.016 

-0.052 

Dy 

0.011 

0.013 

0.041 

0.068 

-0.335 

Ho 

0.001 

0.003 

0.009 

0.004 

-0.075 

Er 

0.003 

-0.028 

0.025 

-0.021 

-0.171 

Tm 

0 

0.001 

0.002 

-0.002 

-0.013 

Yb 

0.001 

0.001 

0.014 

-0.011 

-0.121 

Lu 

0 

0 

0.001 

-0.002 

-0.009 

Eigenvalue 

1292271 

96 

88 

23 

3 

% Var 

100 

0 

0 

0 

0 

Cum % 

100 

100 

100 

100 

100 
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5. Upper background threshold (UBT) 

The upper background threshold (UBT] values are calculated 
for the REEs and the values are presented in Tables 2. The 
U BT is the upper limit of element concentration above which 
the element is considered anomalous. Most of the REEs 
calculated fall below the background values (Table 1, 2] 
except LREEs such as La (samples EB02 to EB015], Ce, Pr 
(samples EB02, EB04 to EB008, EB14 and EB15], Nd 
(samples EB02 to EB09, EB11 to EB15], and Ga (samples 
EB05 to EB07] that fall above the UBT. 

6. Discussion 

6.1. Lanthanide Concentrations 

Lanthanides are lithophile refractory elements and their 
relative proportions in the silicate earth are identical to the 
proportions found in the carbonaceous chondrites [6], 

The studied lanthanides have Ce, La, Nd, Pr enrichment and 
the proportion enrichment of REEs in the study area is still 
remained similar with increasing atomic number. The order 
of the average concentrations of REEs in the study area is 
different to the Oddo-Harkins rule [46], [47] that is Ce> La> 
Nd > Pr > Sm > Gd > Dy > Er > Yb > Eu > Tb > Ho > Tm > Lu. 
Lanthanide concentrations from Bwanebwa area is also 
different to the order of magnitude as that in the earth’s 
crust [8], in the upper continental crust [9], South China Sea 
[48], mean REE values in oceanic crust [49], Redang Island 
Marine Sediment in Terengganu Coastal Waters [50], Bayan 
Obo ores [51] and Edea area [2 7] situated within the Nyong 
Series in Cameroon. These observations show that REEs set 
up a complex group [52] and [27], 

6.2. REE Differentiation Patterns 

Crust normalized REE patterns were used to evaluate REEs 
fractionation. Accordingly, the chondrite-normalized REE 
distribution patterns for the studied stream sediments are 
generally identical, indicating the consistency of geochemical 
distribution of REEs in stream sediments. The curves of the 
REE patterns extend downward from left to right 
characterizing the LREEs enrichment and HREEs depletion. 
The ratio of LaN/YbN quantifies the inclination of the 
normalized REE curves. This ratio of LaN/YbN (11.74-27.93] 
is greater than 1, the curves of LREE incline to right side, 
meaning that the stream sediments is rich in LREEs and low 
in HREEs. The enrichment patterns of LREEs in the studied 
stream sediments is due to the fact that HREEs are more 
soluble and more complex than middle or light REE and are 
more strongly absorbed on most substrates [53] and [27], 
Generally, higher concentrations of LREE are observed in 
soils that developed on phosphate and carbonate rocks, 
whereas the basalt-weathered soils show enrichment in 
HREE [54], The average values, solely for LREEs (616.49 
ppm] and Gd (68.88 ppm] are higher than those reported in 
the UCC: LREEs = 27.86 ppm; Gd = 4.66 ppm [55],These 
average values are also higher than those reported by [56] in 
Brazil. 

6.3. Ce and Eu anomalies 

A very slightly positive Ce anomaly (Ce/Ce* = 1.04-1.07] and 
a strong negative Eu anomaly (Eu/Eu* = 0.097-0.17] were 
observed, indicating that differentiation occurred between 
Ce, Eu and other lanthanides. The depletion or enrichment of 
Ce and Eu usually occurs in the natural environment, which 
may be linked to their oxidation state and mobility under 
different oxidation-reduction conditions [57], Ce of both 
oxidation states are found in soils, but under redox 


conditions Ce 3+ is more easily oxidized to Ce 4+ with higher 
oxygen fugacity and is much less mobile resulting in positive 
Ce anomaly (Ce/Ce* > 1], In addition, the positive Ce 
anomaly recorded in the study area represents an 
enrichment of Ce compared with its neighboring elements. 
Chen et al. (2007] reported that varying degrees of positive 
Ce anomaly was shown in red soils from southern China, 
which reflected the differentiation of Ce caused by different 
oxidation environments during the process of soil-forming. 
Eu is an incompatible element in the trivalent form (Eu 3+ ] in 
an oxidizing magma, but is preferentially incorporated into 
plagioclase in its divalent form (Eu 2+ ] in a reducing magma. 
This ion-exchange process is the basis of the negative Eu 
anomaly (Eu/Eu* < 1], 

6.4. Geochemical Parameters 

The properties of REE composition and differentiation 
patterns are not influenced by weathering, transportation 
process, sedimentation and diagenesis. Meanwhile the 
provenance information carried by REE remains essentially 
unchanged, and consequently REEs could be used as a 
provenance indicator in the geochemical studies [59], The 
binary cross-correlation plots of various geochemical 
parameters were used to discriminate natural variation in 
REE concentrations from other sources which would 
influence REE levels. The cross-plots between total REEs and 
REEs differentiation characteristics, such as LREE/HREE 
ratio, Eu anomaly and Ce anomaly were shown in Fig. 5. The 
geochemical parameters of REE composition from the 
present study exhibit a wide range of variability. The stream 
sediment samples from Bwanebwa area displayed high 
LREE/HREE ratios, distinct negative Eu anomalies and very 
weak positive Ce anomalies when compared with the fluvial 
sediments and REE minerals which influenced by the 
original source of REEs and significantly enhanced by the 
bare tailings [42], Wang et al. (2011] and Xie et al. (2014] 
demonstrated that fine grain size may contribute to the 
enrichment of the REE abundance. There have also been 
reports suggesting that a positive correlation existed 
between the REE differentiation and the gradation of the 
dust grain size [60], These results confirmed that the REE 
distribution characters reflected the compositions of 
materials in the source regions, but are also influenced by 
weathering and pedogenesis. It is worthy to note that soil 
properties such as clay minerals, organic matter, Fe content, 
carbonate, and pH play key roles for the distribution of REEs 
[61]- 

6.5. REE Distribution and Origin 

The concentration of REEs in the study area ranged from 0.1 
to 3758 mg/kg with a mean value of 272.2 mg/kg. Therefore 
the measured mean REEs levels in the stream sediment 
samples were higher than the background values of REEs in 
soil of Baotou region [62], Result observed for the sum of all 
REEs (FREE =1471.56-7674.2 mg/kg, av. = 3811 mg/kg] in 
the present study was also higher than that reported in soil 
(181 mg/kg] by [26], in Japan (98 mg/kg] [63], Australia 
(105 mg/kg] [64] and Germany (305 mg/kg] [65], LREEs 
(La, Ce, Pr, Nd] and Ga fall above the UBT limit and are 
considered as anomalous. A five-factor model of the REEs 
that accounts for 1.998% of the total data cumulative 
variance has been recorded in the study area and the 
significant loading of factors suggests REE fractionation. The 
communality values vary from 0,996 to 1, with a higher 
value reflecting stronger community structure. Factor 1, 
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factor 2 and factor 3 are more important than factor 4 and 
factor 5. The correlation matrix has shown an unusually case 
where REEs correlated higher positively with their 
respective correspondents. These findings confirm the 
aforementioned results and provide evidence for similar 
input sources and common geochemical characteristics of 
the elements. Therefore, lanthanide concentrations in the 
study area were of geogenic origins and not influenced by 
anthropogenic sources. The abundance of REEs in 
Bwanebwa area is related to the lithology made up of 
granites and associated pegmatites. Since lanthanides in 
stream sediments are mobile, they can continuously 
accumulate in stream various pathways such as atmospheric 
deposition [1], mining activities [66] and application of REE 
fertilizers [67], 

Cerium shows the highest average content of REEs, followed 
by La, Nd, Pr, Gd and Snr. This result can be observed in the 
boxplot (Fig. 6) were the median value of Ce is 1792 mg/kg 
with the range of 703 to 3758 mg/kg, La (835 mg/kg; range: 
333-1753 mg/kg], Nd (704 mg/kg; range: 270-1441 mg/kg], 
Pr (191 mg/kg; range: 76.1-392 mg/kg], Gd (66.8 mg/kg; 
range: 27.1-138 mg/kg]. Therefore, patterns of Ce 
distribution show higher than the other lanthanides. The 
REE fractions of these profiles also showed an enrichment of 
LREE relative to the HREE (Fig. 6], These results indicate 
that the REEs patterns of the stream sediment samples 
around Bwanebwa area are consistent with those of Bayan 
Obo ores [51]. Therefore, the results of this work serve as 
guide for lanthanides exploration in Bwanebwa. 



Figure5. Cross-correlation plots of various 
geochemical parameters for REEs 



La Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb lu 


Figure6. Boxplot of rare earth elements from 
Bwanebwa area 

7. CONCLUSION 

Largely distributed around Bwanebwa area, the chemical 
analyses revealed a wide variability of lanthanide 


concentrations with XREE = 1471.56-7674.2 ppm (av. = 
3811 ppm] showing a normal distribution pattern. 
Lanthanides show a high enrichment in La (332-1750 ppm, 
av. = of 863.8 ppm], Ce (703-3758 ppm, av. = 1846.93 ppm], 
Pr (76.11-392 ppm, av. = 196.89 ppm] and Nd (269-1440 
ppm, av. = 715.53 ppm], and the proportion enrichment of 
REEs in the study area is still remained similar with 
increasing atomic number. The order of the average 
concentrations of REEs in stream sediments around 
Bwanebwa area is Ce > La > Nd > Pr > Sm > Gd > Dy > Er > Tb 
> Ho >Yb > Eu > Tm > Lu, which was different to that found 
in the Oddo-Harkins rule (Ce > La > Nd > Pr > Sm > Gd > Dy > 
Er > Yb > Eu > Tb > Ho > Tm > Lu], earth’s crust, upper 
continental crust, South China Sea, mean REE values in 
oceanic crust, Redang Island Marine Sediment in Terengganu 
Coastal Waters, Bayan Obo ores, Edea (Nyong Series], 
suggesting that REEs set up a complex group. 

The LREE/HREE ratio (31.57-35.24, av. = 32.82] showed that 
the content of LREE (1425-7449.8, av. = 3699] is 
significantly higher than that of HREE (46.56-224.4, av. 
112.07], The LREEs account for 97.06% of the total REEs, 
this percentage is in agreement with the trend that was 
observed in the earth’s crust and in the Bayan Obo ores. 
Chondrite-normalized REE distribution patterns are 
generally identical, indicating the consistency of geochemical 
distribution of REEs in stream sediments. Patterns are 
strongly fractionated (LaN/YbN = 11.22-17.92], with LREEs 
significantly fractionated (LaN/YbN of 11.74-27.93] relative 
to HREEs (GdN/YbN = 7.40-16.69, av. = 13.35], implying that 
HREEs are more soluble and more complex than middle or 
light REE and are more strongly absorbed on most 
substrates. Slightly positive Ce anomaly (Ce/Ce* = 1.04-1.07] 
and strong negative Eu anomaly (Eu/Eu* = 0.097-0.17] were 
observed, indicating that differentiation occurred between 
Ce, Eu and other lanthanides. Therefore the positive Ce 
anomaly recorded represents an enrichment of Ce compared 
with its neighboring elements. 

A five-factor model of the REEs that accounts for 1.998% of 
the total data cumulative variance have been observed, with 
the communality values (0,996-1] reflecting stronger 
community structure. PCA results in five components that 
explain 100% of the total variance (eigenvalue], with the 
placement of elements in three major principal components, 
in which elements scattering in the rotated space displaying 
their mutual association. The results provide evidence for 
similar input sources and common geochemical 
characteristics of REEs. Therefore, lanthanide concentrations 
in the study area were of geogenic origins and not influenced 
by anthropogenic sources. The abundance of REEs is related 
to the lithology made up of granites and associated 
pegmatites. Cerium shows the highest average content of 
REEs, followed by La, Nd, Pr, Gd and Sm, hence pattern of Ce 
distribution shows higher than the other lanthanides. LREEs 
(La, Ce, Pr, Nd] and Ga fall above the UBT limit and are 
considered as anomalous. Therefore, the results of this work 
serve as guide for lanthanides exploration in Bwanebwa 
area. 
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